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ABSTRACT 


A first-phase, perturbation analysis is presented for the tracking requirements 
necessary for the measurement of the second-order redshift in a heliocentric probe 
experiment. The clock performance required for such an experiment is also investi- 
gated. 
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THE SECOND-ORDER GRAVITATION REDSHIFT 


- TRACKING ANALYSIS STUDY - 

Jack Jaffe and Robert F. C. Vessot 
Center for Astrophysics 

Harvard College Observatory and Smithsonian Astrophysieal Observatory 
Cambridge, Massachusetts 02138 

1. INTRODUCTION 

The second-order redshift experiment is one of two tests that is currently feasible 
of the second-order term in the theories of gravitation. To date, the only second- 
order test is the measurement of the perihelion advance of a planet or space probe; 
all the other current experimental tests of relativity measure first-order effects 
only. 

A second-order redshift experiment basically consists of sending a space probe 
containing a very stable atomic clock as close to the sun as possible and comparing its 
frequency with an. identical clock on earth by use of light or microwave signals. If the 
gravitational potential difference is sufficiently great, compared to the relative stability 
of the two clocks, the second-order term should be measurable. An analysis of the 
basic theory and equations of this experiment, as well as model orbit calculations, 
have been completed in an earlier phase of this study (Jaffe and Vessot, 1973). 

It is worth noting that one advantage of an experiment of this land is that "idealized” 
orbits, which would necessitate "drag-free" probe tolerances, are not required. 

This paper is basically concerned with a first-phase analysis of the actual tracking 
requirements necessary for measuring the second-order redshift. This analysis has 
then been applied to various model orbits of heliocentric probes, such as the proposed 
NASA— ESRO heliocentric satellite mission. 

This work was supported in part by grant NGR 09-015-205 from the National 
Aeronautics and Space Administration. 
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2. SECOND-ORDER GRAVITATIONAL REDSHIFT 

We can establish the maximum requirements for the necessary accuracy of the 
tracking parameters by application of the perturbation technique. The determination 
of the minimum allowable requirements would necessitate a more extensive analysis, 
by using the correlation and covariance matrix approach. This is beyond the present 
scope of this phase of the study and would be the next logical step in a more detailed 
study of the redshift experiment. 

We will first establish the tracking-accuracy requirements for the parameters as 
they are originally written in the heliocentric coordinate system. These accuracy 
requirements will then be converted into the earth-based coordinate system. 

For an earth-based observer, the "doppler- canceled" redshifted signal z 
received at the earth station from the probe is, in heliocentric coordinates (Jaffe and 
Vessot, 1973), 



(1) 


2 



9 2 2 2 2 
where m = GM/c , p. = v./c (G and c are set = 1 throughout), p. = p + r"p 0 , 

11 i i 

A sp! • 77 0 , and A' = (f 0 • T , in which T, = [/A/B Vl - (Bi 2 /r 2 )> is the 

G 1 p £» . 

generalized "direction vector" of the light signal from position i to position j; i is a 

2 2 

constant of the motion for the light signal, B = 1 - (2m/r) + (p - -y) (2m /r ), and 
A = 1 + -y (2m/r). The unsubs cripted parameters, p and y, are the so-called 
Eddington-Robertson coefficients (Jaffe and Vessot, 1973). 



3. DETERMINATION OF 1 IN TERMS OF THE BOUNDARY CONDITIONS 


Before proceeding with the specific perturbation analysis, we must first determine 
i in terms of r., 0 1} r , and 0 . This can be done by using the geodesic equations of 
motion for light. We first carry out the calculation in the absence of gravity, and then 
we repeat this analysis for the spherically symmetric gravitational field. 

In the absence of gravity, all the g„ = 1, and the relevant geodesic equations for 
light in the (r, 0, tj>, t) coordinate system are as follows (Jaffe, 1969): 

2 d9 , , _ 

r ds = constant = P > 

i = constant =q , 

2- = constant ; 

q 

i can be interpreted as the impact parameter as measured at infinity.' The <f> equation 
is not necessary here, for light can be confined to motion on a plane; this is true both 
in the classical and in the general relativistic cases. 

For g.j = 1, the metric is simply , 

, 2 ,.2 , 2 2 ,„2 
ds = dt - dr - r d9 . 

2 

For the case of light, ds = 0. 

Combining these equations, we have 
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The ± sign depends on the direction of the light signal. This equation can be integrated 

r , • f , < ’ f ■ 

to yield 


9 = ± cos 1 + 9 0 , 


dr 

where 0 Q is the point at which 


= 0. The corresponding value of r, the so- 


9=9 


0 


dr 


called turning point, is defined as d^, i.e., ^ 


= 0 . 


r=d^ 


We have 


9 1 " cos 1 (r 1 ) + 9 0 » 


°2 = C0S_1 ( 4 ,) + 9 0 ' 


Subtracting the second equation from the first and taking the cosine of both sides yield 

„ Q v 1 r i)2 .2 ,2.1/2 , 2 .2,1/2, 

cos (0 2 - 9 X ) = [ * ~ (r J _ - i ) ' (r 2 - i ) ' ] . 
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By solving for f , we finally have 


i = 



r l r 2 Sin < 9 2 - 9 1 ) 

+ r^ - 2r^r 2 cos {0 2 



( 2 ) 


Since, classically, light moves in a straight line, the "impact parameter at 
infinity, " i, is indeed the value of the distance of closest approach at the point of 
closest approach. Thus, 
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This is not the case in general relativity, as will be seen below. 


We now repeat the calculations in the general relativistic framework. We will 
use the exact Schwarzschild field here, instead of the Eddington-Robertson expansion, 
for the spherically symmetric gravitational field, because it is somewhat simpler to 
manipulate and is quite sufficient for our present purposes. The Schwarzschild field 
is given as 


ds 2 = A dt 2 



2 JQ 2 

r d0 




where A=l - 2m/r. For light, ds Z = 0. 


The relevant geodesic equations are 

2 d6 , . _ 

r ^ = constant = p , 

A — = constant = q , 




where f can again be interpreted as the impact parameter as measured at infinity. 
Combining these equations gives 


, . \ 2 4 

/dr \ _ r 

\ d 9 / “ .2 


A r 


( 3 ) 


t* 

^0 = 0 at r = d. This yields 


i 2 =* 


A 1 - (2m/ d) 



or 


yi - (2m/d) 


( 4 ) 


Thus, i does not have the same value as the turning point, d. By integrating 
equation (3), we have 



- d 3 /2 f 


d»m 




dr 


~ w 3 . ,3 _ . 3 ,3, 

r ' fr d-d r - 2m (r - d ) 


= d / 


dr 


■v; 


*. 2 ,2 L 

r * r - d 


1 + 2. + md 


d r(r + d) 




- COS 


■‘H V r 2 - d2 [^ri] + % > 


dr 

where is the value of 9 where -jg 


= 0 . 


9 = 9 , 


Accordingly, 

n _ -1 d , m ./ 2 71 

0 A - 8 = cos — + — w r - d 
0 1 r, d t 1 


2r^ 4- d 
r 1 (r 1 + d) 


-Id m 

J 2 r ,2. 

" 2r 2 + d 1 

COS + "T“ 

r 2 d 

V r 2 d 

_ r 2 ( r 2 +d)_ 


(5) 
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Adding these two equations and taking the cosine of both sides yield 


/ -1 d . 

cos ( cos — + cos 

V r 2 


* ft 17 


cos (9 2 - e x ) + ^ 




2 2 2r 2 +d 0 

r‘-d? 2 — 


2 0 r 2< r 2 + d 0 > 


JT ,2 2r i + d o 

» r i - d 0 


l< r l + d o> 


sin (0 2 - 0 : ) 


Solving for d, we obtain 


1/2 


r 2 r 1 |l - £cos (0 2 - 0^ + (m/d 0 ) E sin (0 2 - 6j)]“j 

— — - ‘ -J/2 

jr 2 + r 2 - 2r 1 r 2 [cos (0 2 - Gj) + (m/dg) E sin (© 2 - 0^ j j 


r 2 r l sin < 9 2 ~ 9 1 ) 




r i + r 2 _ 2r i r 2 COS (0 2 ' °1 ) 


1 - ^ E 
d 0 


ctn (0 2 - 9j) 


r l r 2 Sin (0 2 ' 9 1 ) 


r l + r 2 " 2r i r 2 COS (0 2 " 9 1 ) -1 


( 6 ) 


where 


E-V7T7 

’ 2 0 r 2< r 2 + d 0> ’ 1 0 r l (r l 


2r , + d. 


+ d 0> 


8 



Comparing equation (6) with equation (2), we see that 
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4. PERTURBATION ANALYSIS IN THE HELIOCENTRIC SYSTEM 

We will first carry out a perturbation analysis with respect to the heliocentric 

space-probe parameters r , 0 ? v , - Errors due to uncertainties in the earth- 

A r<2 ^2 

based observer’s parameters (r^ 0p . . . ) will be significantly smaller, and need not 
be considered here in an analysis of the most stringent tracking accuracies required. 

4. 1 Radial Velocity 

Differentiating equation {1) with respect to 6 , we have 

2 



Lower order terms in v/c need not be retained here, since we are solving for the 
maximum uncertainty in the parameter that will still allow detection of the second- 
order (in m/r = fourth order in v/c) redshift effect. 

4.2 Angular Velocity 

Differentiating equation (1) with respect to (3_ , we have 

°2 

Az = <- r^ (3 0 + A g l ) Ap 0 . (10) 

2 2 

4. 3 Radial Distance 


We differentiate equation (1) with respect to r^. This result is more involved than 
the preceding equations, since i is also a function of r^ and 0^. We find 
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Az 





Ar n + A 
2 e 






where 


£ 8 



r l r 2 sin (9 2 ~ 9 1> 

+ v\ - 2r x r 2 cos (9 2 - 0j) 


( 12 ) 


(The corrections in m/d to £ need not be retained here, for the same reasons as given 
above. ) And, 

= I l Sin(e 2~ 9 l )lT l- r 2 C ° B(e 2- 9 l >l .... 

3r 2 [r 3 ~ t 2 2 - cos (0 2 - e~iT 3/2 

4. 4 Angle 


Differentiating equation (1) with respect to 0 2 yields 



( 14 ) 
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where 


»1 r i r 2 C0S < 9 2 - e i> [r l + r 2 - r l r 2 C0S < e 2 " 9 1 )] - r l r 2 

90 2 [ Tj + r 2 " 2r i r 2 cos ( 0 2 " e i^ 3//2 


(15) 
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5. CONVERSION FROM HELIOCENTRIC SYSTEM TO 
GEOCENTRIC SYSTEM 

It is sufficient at this stage of the analysis to make a plane -geometry conversion 
from the sun-centered coordinate system to the earth-centered system. From 
Figure 1, we have 

r~= R-t-f , 

r 2 = R 2 + p 2 + 2R • fT 


and 


f - R+ p , 


j3 2 = f 2 = R 2 + p 2 + 2R * p . 
c=l 


The unit vectors in the r and 0 directions can be written as 


e r = cos 0 i + sin 9 j 


= Vl- sin 2 0 i + sin 0 j , 


0 ~ - sin 0i.+ cos 0 f 


- - sin 0 


l + ^1 - sin 2 0 j , 


where l and ^ are the unit vectors in the x and y directions (in the heliocentric system). 
By the law of sines, 


sin 9 sin <t> 


( 16 ) 
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or 


sin 0 - £ sin 6 „ 

r Y 

We finally have 

e r = ^1 - sin 2 4> i + £ sin 4> j , 

with 

- 2Rp cos <f> , 

, A 

an expression for written solely in terms of the earth-based coordinates p, <f> } and 
R. Similarly, 

= " r sin ^ * + V 1 ~ % sin2 $ J • 

r 

We can now write P y and p 0 completely in terms of the earth-based coordinates: 
P r ~ P"’ e r = (R+ p) ♦ e r , 
where p has the components 

P = P P + P <(> f 

= - (p sin <f> + p cos <f>) i + (p cos <j> <p - p sin <j>) j . 

The earth T s motion, R, is assumed here to be uniform and circular; i. e., 

R = - R to sin cot i + R co cos cot f . 

Corrections to this assumption will be discussed in a later section. 
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We finally obtain 


P r = (R+ P) • e r 


2 

= - (R to sin cot + p sin <(> <j> + p cos (}>) V 1 " o s * n ^ 'f 5 


r. 


+ (R co cos cot + p cos <f) <£ - p sin <j>) £ sin <j> , 


( 17 ) 


with 


2 2 2 

r = R + p - 2Rp cos 4> • 


Similarly, 


% = (R+p) • e 0 


— (R <o sin cot + p sin tj> 4> + p cos <£>) £ sin <j> 


+ (R to cos tot + p cos 4> $ - p sin <f>) ^ 4 1 


2 o 
p . 2 

— am 


(18) 


(19) 
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6. PERTURBATION ANALYSIS IN THE GEOCENTRIC SYSTEM 


The preceding section has established the functional relationships of the param- 
eters in the heliocentric system with respect to those in the geocentric system. We 
must now consider the resultant uncertainty in these geocentric parameters that is 
associated with the allowable uncertainty in the heliocentric parameters. 


From equation (18), 


r 2 = (R 2 + p 2 - 2Rp cos <f>) 1//2 


Therefore, 


9r ? 9r ? 9r 2 

* R+ 5r Ap + sr A|> 


R - P cos <fr 

<R 2 + p 2 - 2Rp cos 


AR 


2 P-Hcos<|> 

(R^ + p - 2Rp cos <j>) v ^ 


R P sin( fr ^ M 

,2 . 2 “ ,, 1/2 


(R + p - 2Rp cos <(>)■ 


(20a) 


(20b) 


(20c) 



From equation (16), 
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Therefore, 


39 90 99 9 


R - P cos <t> P sin <f> 

fl t 2,\ 2 " 1,172 3^^ 

l 1 “ (P /r 9 ) sin 4>] ' r 


[1 - (pW) sin 2 4>] 1//2 


sin cb p sin d> . „ „ „ ,, 

~zr~*~ ~ 3 r (p - R cos <f>) 


L 2 r 2 


Ap 


[1 - (p 2 /r 2 ) sin 2 cf>] J / 2 


P c - s ^ - ( Rp S i n <|,) 

r 2 r^ 


Ac|> 


P 2 

r r 

2 


From equation (17), 


3P r 3P r 9P r 3P r; 8 P r 9 P r 
^r = 3^ Ap + 7T J ' A P + 9R 1 AR+ d^ A 

2 r C7 P C7<p 


. I ~2 ~ 

sin 4> cj> Jl - g sin 4> 

1 


* 

- (R co sin cot + p sin <}> <|> + p cos t)>) 


X | 1 - -^C sin 2 4> 


r 1 / 2 r 


2 2 2 
P sin <i> . p sin c|> , „ _ ±v 

- z — — ii + X — x (p - R COS 4>) 

r r 

2 2 


+ cos 4> (£>-£• sin <J) + (R co cos cot + p cos 4* — P sin 4>) 

r 2 


X 


'sjai _ (p _ R 00s # 


Ap 


(21a) 


(21b) 


(21c) 


(22a) 
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cos 4) 


^ 1 - ^ sin 2 ^ + sin 2 4> 


(22b) 


- ( p cos 4> <f> - p sin 4>) 


yl-^2 sin2 4> “ ( R ^ sin wt + P sin $ ^ + P cos $ 


X 


9 3 2 

2p sin 6 cos <|» 2p sin 4 R 

2 4 

r 2 r 2 


- (p sin <f> + p cos <j)) ^ s ^ n ^ 

T 


+ (R o) cos cot H" p cos <() <j> - p sin <f>) 


X r p sin <{> J | Af> 


(22c) 


p sin <f> 


F 


sin^ <j> 

A 

r 2 


p cos cf) sin cj> | A4> 
r 2 


(22d) 




l . 2 

— enn 


co sin cot V 1 “ 2 sin <j> 
rt 


o 2 

+ (R w sin oot + p sin <)> <j> + p cos 40 1 1 “ jj sin ^ 

r 2 


-1/2 


2 2 

X P s ^ n ^ (R - P cos 4 ) - (0 cos tot sin 4 


+ (R to cos cot + p cos 4 <j) - 0 sin 4>) (R - p cos 4) 


AR (22e) 
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(R sin cot + co R cos cot) 


, P_ . 2 , 
1 - o sin <j> 


9 P 

+ (R cos cot - oo" R sin cot) — sin cf) Aco 


From equation (19), 


9 \ *e 2 \ 9 \ . 8p e 2 8p 9 2 

^e 2 = 8T 1 * + AR + — A “ 


- < sin <f> cj) *j"~ sin 6 + (R co sin cot + p sin 4> <j> +. p cos <£) 


sin 


-g — (p — R COS cj>) 4- COS <4> cj> 


1 P_ . 2 , 
1 - 1 ~2 sin <j> 

r 2 


+ (R oo cos cot + p cos <f> cj> - p sin cf>) 


! i i+£ i^£i (p . Rcos4 , ) 


— (R co sin cot + p sin 4> + p cos <j>) sin <jo 


i £_ • 2 

1 - x — sm 


+ (R oo cos cot + p cos <f> <j> - p sin cj>) 


X 2 y RcOS + ^ Ap 

(R Z + - 2Rp cos <j)) 1/4 \ 
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(23b) 


f- 


2 9 

+ — ( cos <j> -£• sin cj> - sin cj> J 1 - ^ sin" <j> ) Ap 


+ — < (p cos 4> - p sin <()) -2- sin <j> 

r 2 / r 2 


• . ..Ip cos <b p sin 4 n , 

+ (R co sin cot + p sin 4> cj> + p cos 43 ) - — 31 - £ — g — 31 R P sin <j> 

\ 2 r 2 


f- 


2 

- (p sin <}> cj) + p cos <b) +1 1 - ^ sin 4> + (R co cos cot + p cos <j> <j> 


X 


2 2 3 

2p sin <fr cos 4> 2p sin <j> R p 

2 4 

y r 

2 2 


(R co sin cot + p sin <{> cj> + p cos <j>) sin cj) 


2 2 

+ (R co cos cot + p cos cj) <j> - p sin cj)) W 1 - ^ sin <j> 


F 


vy R p sin cj> v a 1 

X o 2 • 172" } A* 

(R -s- p - 2Rp cos <j>) ' 


+ ~ J P sin^ 4 ^7- + p cos cj) y 1 - sin^ <j> ] Aj> 
2 \ 2 r 2 


p sin cj>) 


(23c) 


(23d) 
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+ — < to sin wt -£• sin d> - (R w sin wt + p sin 4> + |5 cos 4>) ^ 5 ) p * (R - p 

r 2 l *2 r 2 


+ to COS cot 


I 2 2 

W 1 - ^ sin*" 1 <}) + (R w cos tot + p cos 4> <£ - p sin <j>) 


2 \ 2 . 2 
1 - sin 2 <j> ) X ^ ^ ^ {R - p cos <f>) 


.L 


(R to sin tot + p sin 4> 4> + P cos tj>) -*— sin <Jj 

r 2 




2 2 

+ (R to cos tot + p cos cj> tj> - p sin tj>) V 1 ” 2 S ™’ J * 


X 


R - p cos (f> 

(R 2 + p 2 - 2R p cos <j>) 


AR 


(23e) 


(R sin tot + to 2 R cos wt) sin cj> 

r 2 


+ (R cos cot - to R sin tot) 


rj_ - 2 , 

y 1 - - J ~2 sin cf> 


Ato 


(23f) 


cos 4>) 
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7. CORRECTIONS TO IDEALIZED MODELS 


The redshift equation, equation (1), and the various other forms derived in the 
preceding study (Jaffe and Vessot, 1973) are completely general; they have been set 
up without any restrictions on (L. Our idealized model here has assumed a "point" 
earth moving in a perfect circle around the sun; the velocity of revolution is 
~30 km/sec. Of course, the earth is not a point, and it is necessary to investigate 
what corrections to our present analysis would be introduced by considering a more 
realistic model of the earth. 

The earth has a mean radius of ~6400 km and rotates on its axis with a velocity 
of rotation of ~0. 47 km/sec. Since this velocity of rotation is about two orders of 
magnitude less than the velocity of revolution around the sun, it is not necessary to 
include the secondary velocity effects in this present determination of the greatest 
permissible uncertainty in the tracking parameters. Similarly, since the earth's 
radius is approximately four orders of magnitude less than the earth— sun distance, 
inclusion of a finite-radius earth in the transformation to a geodetic system would 
negligibly change the resultant equations here. 

All secondary effects, such as the earth’s rotation, motion about the earth— moon 
barycenter, and the earth and moon's actual gravitational potential, can be included 
in the final working program, by use of the completely general formalism already 
derived. As was emphasized in the preceding study (Jaffe and Vessot, 1973), all 
coordinates, parameters, and constants must eventually be written in an operational 
framework, completely in terms of observables (invariants), in order to eliminate 
any ambiguity in the interpretation of a coordinate system. 
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8. ANALYSIS OF TRACKING-ACCURACY REQUIREMENTS FOR 
HELIOCENTRIC PROBE MODEL ORBITS 

The perturbation equations derived in the preceding sections were programed 
and applied to various model orbits of a heliocentric probe, such as the proposed 
NASA— ESRO heliocentric satellite mission. 

The general analysis of the tracking- accuracy requirements is as follows: For a 
given time in the orbit, all the various -order terms in the redshift were calculated. 

The smallest term was then used to determine the maximum allowable uncertainties 
in the heliocentric parameters, as given by the equations in Section 4. These uncer- 
tainties in the heliocentric parameters were then converted to corresponding uncer- 
tainties in the geodetic parameters, through the equations in Sections 5 and 6. 
Obviously, an uncertainty in a geodetic parameter such as p could correspond to 
uncertainties in more than one heliocentric parameter, in this case, r 0 , 0 O , P , ; 

“ “ r 2 ”2 

all these possibilities were computed, commensurate with the actual allowable 
uncertainties in the heliocentric parameters, and the program then chose the maximum 
uncertainty in each case. The results are shown in Figures 2 to 4. 
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9. STATE OF THE ART IN DEEP-SPACE-PROBE TRACKING 


The present state of the art in deep-space-probe tracking is illustrated in Table 1, 
which gives the 1-sigma accuracies achievable with observations carried out both with 
present-day and with projected future capabilities in instrumentation and data-process- 
ing methods. The given values do not take into account the possible contributions to 
the achievable maximum accuracy coming from a priori knowledge of the space 
trajectory. 

All figures apply to dual -frequency observations, which remove the effects of the 

interplanetary plasma as well as those of the earth’s ionosphere. The angular accuracy 

values apply to earth-based very long-baseline interferometry (VLBI) systems, in 

which tropospheric effects have been removed (down to 0. 1%) by columnar radiometric 

observations (conducted at water vapor and oxygen absorption lines) of the tropospheric 

path along the direction of, and simultaneously with, the VLBI beam. The angular 

—9 —10 

accuracy figures could be substantially lowered (to 10 or even 10 rad) by using a 
100, 000 -km VLBI baseline established between two satellites, thereby totally removing 
the effect of the time-variable and space-variable earth’s atmosphere. With this 

• _i3 

spacebome VLBI system, the angular rate accuracies could be improved to 10 or 
10 rad/sec. 

Table 1. Achievable Accuracies with State -of-the -Art Deep-Space -Probe 
Tracking Techniques 


Capability 

Radial 

range 

<m) 

Radial 

velocity 

(mm/sec) 

Angle 

(rad) 

Angular 
velocity 
(rad/ sec) 

Present 

16 

0.5 

2 X 10" 7 

io - 12 

Future 

1 

0.05 

2 X 10~ 8 

io - 13 
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10. COMPARISON OF REQUIRED AND CURRENTLY 
ACHIEVABLE TRACKING ACCURACY 


The analysis indicates that the current tracking capability is sufficient to determine 
the range p and earth-sun distance R accurately enough for a second-order redshift 
experiment. The JPL Deep Space Network (DSN) tracking program, for example^ has 
a present range tracking accuracy of approximately 16 m and a projected, future 
accuracy of about 1 m. 

With regard to range rate p, for all model orbits considered, the present con- 
ventional tracking schemes are about two orders of magnitude too crude for some parts 
of the probe orbit, but nearly sufficient near perihelion (which corresponds to the 
greatest possible redshift with respect to the earth). The projected, future accuracy 
is still a bit less than is required over the entire orbit. 

The angular tracking accuracy is also too crude in the regions far from perihelion 
by roughly three orders of magnitude. Proposed accuracy for earth-based terminals 
is still roughly two orders of magnitude too crude in these regions, even with the 
assumption of the best estimate for earth-based VLBI tracking in the near future. As 
mentioned in Section 9, a VLBI tracking system using satellites in earth orbit, although 
a somewhat distant possibility, would be one way, for example, to achieve the required 
tracking accuracy. 

The angular velocity and u requirements are quite stringent for all model orbits. 
This question requires more detailed study and depends on the specific tracking 
scheme to be used. 


25 



11. CLOCK ACCURACY REQUIREMENTS 


We can. determine the necessary performance required of the clocks for a test of 

4 

the second-order redshift term by plotting the minimum (3 term in the redshift as a 

function of the orbital time. Figures 5 to 7 indicate that the requirements are least 

stringent near perihelion, where the redshift is maximal. (The curves are, of course, 

not necessarily symmetric around perihelion, since the relative juxtaposition of the 

probe and earth at different times will cause the |3 or m|3 terms to become as 

2 

significant as the pure m term.) 

The present state of the art in atomic clocks is advancing toward stability in the 

“16 5 g 

10 region, for averaging intervals of the order of 10 to 10 sec. Ground-based 

15 

maser devices have already shown stability of 2 parts in 10 for time intervals of 10 to 
4 

10 sec. Development is in progress of a space-qualified hydrogen-maser clock for 
a terrestrial probe experiment; development of space -qualified cesium and rubidium 
devices is also expected in the next few years. 
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12. RECOMMENDATIONS 

This perturbation analysis is intended to present the "worst case" situation. It 
has maximized the tracking requirements for the second-order redshift experiment. 
Correlations of the parameters with respect to the earth and probe orbits, which 
would probably tend to make the actual requirements somewhat less stringent, have 
still to be considered. This further study requires more extensive effort and is 
recommended as the next step in the analysis of a second-order redshift experiment. 
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SPACECRAFT 


Figure 1. Tracking geometry. 
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Figures 2, 3, and 4. Allowable uncertainties in the geodetic parameters, as a function 
of the probe period. A plot of the earth— sun distance uncertainty AR, and the earth's 
angular velocity of revolution Au, have been omitted for simplicity; they are either of 
the same order, or approximately 50% more stringent, than Ap and Ap, respectively. 
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Figure 3. 
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Figures 5, 6, and 7. Minimum term of order c , as given by equation (1), as a 
function of probe period. 
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Figure 5. 
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